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Abstract 

We enumerate self-avoiding walks and polygons, counted by perimeter, on the 
quasiperiodic rhombic Penrose and Ammann-Beenker tilings, thereby considerably 
extending previous results. In contrast to similar problems on regular lattices, these 
numbers depend on the chosen initial vertex. We compare different ways of counting 
and demonstrate that suitable averaging improves convergence to the asymptotic 
regime. This leads to improved estimates for critical points and exponents, which 
support the conjecture that self-avoiding walks on quasiperiodic tilings belong to the 
same universality class as self-avoiding walks on the square lattice. For polygons, the 
obtained enumeration data does not allow us to draw decisive conclusions about the 
exponent. 

1 Introduction 

Quasiperiodic tilings are most widely used for the description of quasicrystals. With ap- 
propriate atomic decorations of the vertices, they serve as structure models which explain 
physical properties of quasicrystals [10]. From a theoretical point of view, they are ideal- 
isations of real substances on which the usual models of statistical physics like the Ising 
model may be studied |2H1 1201 12^1 • Quasiperiodic tilings arose before the discovery of 
quasicrystals, however, more as an object of aesthetic interest in geometry fI7\ l2Uj. 

From a combinatorial point of view, they provide an interesting example of a non- 
periodic yet structured graph where typical problems of combinatorics like the counting of 
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objects on a lattice become more complex. This is fundamentally different from counting 
problems on semiregular lattices, where the underlying translational invariance is still 
present flTl, and also from self-similar graphs, where the self-similarity allows for the 
solvability of some counting problems. 

Consider for example the counting problem of n-step walks on a quasiperiodic tiling. 
This number depends on the chosen starting point. On a lattice, this phenomenon does 
not occur due to translation invariance. Questions arise, such as if the universal properties 
of the walks like critical exponents j^Hl are changed for quasiperiodic tilings^ , and how 
different ways of counting affect asymptotic properties. The first question has been inves- 
tigated for self-avoiding walks (SAWs) and self-avoiding polygons (SAPs). These are walks 
or loops which do not visit the same vertex twice. Extrapolation of exact enumeration 
data for a number of quasiperiodic tilings |2H1 indicates that the critical exponents 7 
for SAWs and a for SAPs are consistent with the corresponding values on regular lattices 
a = 1/2 and 7 = 43/32. In [22j, the related problem of SAWs on Penrose random tilings 
[T^ has been studied by Monte Carlo simulations, the results suggesting the same mean 
square displacement exponent as in the lattice situation. The studies |2H] suffered how- 
ever from strong finite-size effects due to the relatively short series data available. This is 
mainly due to the fact that the finite- lattice method jH] , being the most successful method 
known to date for walk enumeration on regular lattices [HEj, cannot easily be applied here, 
and so we used the slower method of direct counting. The problem in applying the FLM 
is discussed in Section 3 below, see also [211 • Another reason for the pronounced finite-size 
behaviour is that the number of walks or loops depends on the chosen starting point. One 
might suspect that suitable averaging over different starting points reduces these effects, 
leading to behaviour comparable to the square lattice case. In this paper, we analyse three 
different methods of counting in detail. Whereas the first one depends on a chosen vertex 
of the tiling, the last two are averages over the whole tiling. 

• Fixed origin walks. We count the number of n-step self-avoiding walks emanating 
from a given vertex. This number depends on the chosen vertex. 

• Mean number of walks. We count all translationally inequivalent n-step self-avoiding 
walks which may occur anywhere in the infinite tiling, weighted by their occurrence 
probability. For tilings with quasicrystallographic fc-fold symmetries, these proba- 
bilities are numbers in the underlying module U\e^'^'^l^\. This leads to a generating 
function which has non-integer coefficients. 

• Total number of walks. Here we count the number of translationally inequivalent 
n-step self-avoiding walks which may occur anywhere in the tiling. This number is 
bigger than the number of fixed origin walks, and by definition, takes into account 
vertices over the whole tiling. 

Note that two self-avoiding walks (polygons) are translationally equivalent iff they have, up 
to a translation, the same vertex coordinates. For self-avoiding polygons, we will employ 

^For ferromagnetic spin systems on quasiperiodic graphs, a heuristic criterion determines whether its 
behaviour is different from the lattice situation |25| . 
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the second and the third method of counting. We do not distinguish between different 
filhngs of the interior of the polygon. For SAPs, the second method has been implemented 
previously to obtain the high temperature expansion of the Ising model. The mean 
number of SAPs up to length 2n = 18 has been determined on the Ammann-Beenker tiling 
IT^ and on the rhombic Penrose tiling IHj. 

We counted SAWs and SAPs on the Ammann-Beenker tiling and the rhombic Penrose 
tiling and compared different counting schemes, thereby extending and generalising the 
previous approaches to counting SAWs [Sj and SAPs [2H]- Generally speaking, averag- 
ing reduces oscillation of data due to finite size effects, providing improved estimates for 
critical points and critical exponents. Within numerical accuracy, we cannot rule out the 
universality hypothesis that SAWs on the Ammann-Beenker and on the rhombic Penrose 
tilings have the same exponents as on the square lattice. The data for the total number of 
walks (polygons) gives a different exponent, reflecting the fact that the number of patches 
grows quadratically with the patch size, in contrast to the lattice case j^Ill^. The limited 
quantity of SAP data does not allow us to draw decisive conclusions about exponents. 

This paper is organised as follows. The next chapters describe the algorithms used for 
the generation of the tilings and for the computation of the numbers of walks and their 
mean values. The following chapter is devoted to the asymptotic analysis of the series and 
to a comparison of the different approaches. This is concluded by a discussion of possible 
future work. 

2 Graph generation 

Quasiperiodic tilings in may be obtained by projecting certain subsets of lattices from 
a higher- dimensional space into W^. This is described by a cut-and-project scheme, 
summarised in the following diagram. 

TTll 7r_L 

E|| ~ M'^' < — E = W — > E^c^ W 

u \ u / u 

1-1 dense 
L\\ = 7r||(L) L lattice W polytope 

It consists of a Euclidean vector space E, together with orthogonal projections tt\\ and 
7r_L. The vector spaces E\\ = vr||(£') and E± = 7i±{E) are called direct and internal space, 
respectively. Let L C -E be a lattice. The projections are such that 71^1^ is one-to-one and 
7!'±{L) is dense in E± (or dense in some subspace of E±). Let W C E-^ be a polytope (or 
a finite union of polytopes). The set W is also called the acceptance window. The set of 
tiling vertices A(l^) is defined by 

A{W) = e L\\\xe L andx^ eW}. (1) 

The edges of the tiling are defined by the following rule: The tiling vertices vr||(x) and 
7r||(y) are adjacent iff the lattice vectors x and y are adjacent. 
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For the Ammann-Beenker tiling ^ ^] , we have n = 4 and d = m 
L = Z^. The projections Trii and tt± are defined as follows. For x G M", 
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The acceptance window W C M™' is a regular octagon with unit side length centred at the 
origin, having edges perpendicular to the axes. A typical patch is shown in Figure HJ 




Figure 1: A patch of the Ammann-Beenker tiling. 



For the rhombic Penrose tiling j27| 16]. we have n = 5, d = 2 and m = 4. The lattice is 
L = 11'. The projections ti\\ and 7rj_ are, for x G M", defined by 
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The acceptance window W C is made up of four regular pentagons in the planes 
Xx3 = 0, 1, 2, 3. The pentagons in the and 3 Xs-planes have unit side length and the 
others have side length 2cos|. Each pentagon is centred at x_|_i = 0, x_|_2 = 0. Pentagons 
and 2 have an edge crossing the positive x_|_i axis at right angles while pentagons 1 and 
3 are rotated through |. A typical patch is shown in Figure El 
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Figure 2: A patch of the rhombic Penrose tihng. 



We remark in passing that a more natural embedding of the rhombic Penrose tihng is 
the root lattice ^4, see also jS], but is more convenient for computations. Moreover, 
the Ammann-Beenker tiling and the rhombic Penrose tiling may be alternatively defined 
by inflation rules for their prototiles |31 ■ 

3 Enumeration 

A self-avoiding walk on a graph is a path, beginning at an origin vertex, which never 
visits a vertex more than once. The SAW on the lattice Z"^ is a well-studied object (see, 
for example, [211 1231 )• The number c„ of translationally inequivalent n-step walks on 
a regular lattice is clearly independent of the choice of origin vertex. Hence this series 
is representative of the entire lattice. The enumeration of SAWs on non-periodic tilings 
introduces complications to the interpretation of the series C{x) = X]n>o'^"-^"- '^^^^ 
because the possible origin vertices produce an infinite range of different series C (x) . Each 
origin produces a different series which is representative only of that vertex's immediate 
neighbourhood in the tiling. The question then is, how do we obtain a SAW series which 
is representative of the whole tiling? In this paper we adopt three different approaches to 
enumerating SAWs on quasiperiodic tilings, as described in the introduction. They are: 

• Fixed origin walks. 

• Mean number of walks. 

• Total number of walks. 
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3.1 Fixed origin walks 

We take a random selection of origin vertices x G L (if x_l ^ the vertex is not a 
suitable choice and is ignored). For each suitable origin, we generate the neighbourhood 
of the vertex, including all vertices up to some Euclidean distance away. Two such 
neighbourhoods are shown in Figure and Figure El We enumerate all SAWs from the 
origin up to length n in the neighbourhood using backtracking This takes time 

proportional to the number of walks c„. Unfortunately, the transfer matrix approaches 
used to enumerate SAWs on regular two-dimensional lattices in less time cannot be used 
on this problem without major adaptation: Since a SAW of A^ steps may reach a vertex A^ 
steps from the origin, we would need to consider every possible tiling patch of radius A^. 
The finite-lattice method's transfer matrix stage would then need to be adapted to each 
tiling patch, or generalised to handle them all. 




x± = (0,0) XX = (1,0) XX = (1-^,^) xx = (l-v^,0) 

Figure 3: The actual Ammann-Beenker neighbourhoods chosen for the enumerations. 




xx = (0,0,1) xx = (0,0,0) xx = (0.5,0.5, 1) x^ = (0.25, 0.5, 0) 



Figure 4: The actual rhombic Penrose neighbourhoods chosen for the enumerations. 

If we apply this method to counting walks on a regular lattice it is clear that we 
would always produce the usual SAW series for that lattice. If each of the series in Table 
HI and Table |21 showed lattice consistent properties, it would be a good indication that 
these properties belong to the entire tiling. The actual neighbourhoods chosen for the 
enumerations are shown in Figure El and Figure |31 for the Ammann-Beenker and Penrose 
tilings respectively. 
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n 


x± = (0,0) 


x± = (l,0) 


xx = (l-^,^) 


x^ = (l-^/2,0) 





1 


1 


1 


1 


1 


8 


3 


4 


5 


2 


16 


13 


12 


16 


3 


48 


34 


46 


42 


4 


144 


108 


108 


152 


5 


448 


292 


374 


388 


6 


1088 


952 


976 


1194 


7 


3680 


2458 


3042 


3412 


8 


9584 


7746 


8330 


9678 


9 


28336 


21348 


24556 


27218 


10 


82960 


61478 


68376 


79150 


11 


225408 


177230 


197820 


217562 


12 


657536 


495808 


554108 


628996 


13 


1834768 


1412152 


1576464 


1741464 


14 


5140752 


3985706 


4400920 


4968606 


15 


14584112 


11125408 


12531794 


13724682 


16 


40222672 


31617786 


34541864 


39209054 


17 


114683280 


87149372 


98846548 


107503768 


18 


313146848 


248799302 


270221012 


306845714 


19 


896810944 


680172768 


773046904 


840463852 


20 


2437468000 


1943692238 


2109562128 


2386875508 


21 


6958267152 


5303535884 


6011045200 


6548653714 


22 


18981078176 


15086983820 


16431248782 


18500898140 


23 


53728620912 


41295324398 


46538635588 


50883461478 


24 


147472084608 


116624466842 


127704810544 


142927122532 


25 


413887940176 









Table 1: The number of n-step fixed origin SAWs for various starting points x_l in the 
Ammann-Beenker tiling, with starting point coordinates {x, y) given in the internal space. 
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n 


x± = (0,0, 1) 


x_L = (0, 0, 0) 


xj_ = (0.5, 0.5, 1) 


xj_ — (0.25, 0.5, 0) 





1 


1 


1 


1 


1 


5 


5 


5 


4 


2 


20 


10 


1 A 

14 


12 


3 


40 


40 


50 


46 


4 


160 


130 


130 


112 


5 


450 


310 


406 


394 


6 


1170 


1140 


1177 


938 


7 


4000 


2680 


3316 


3416 


8 


9480 


9360 


9723 


7866 


9 


32910 


23150 


27356 


27312 


iU 


( Duyu 


i zozu 


777/1 7 


uOiOU 


11 


262250 


196980 


224102 


215924 


12 


619460 


555290 


615549 


545062 


13 


2050500 


1646990 


1812802 


1698548 


14 


5052310 


4292010 


4869786 


4411293 


15 


15828550 


13403280 


14455725 


13367278 


16 


41103090 


33637420 


38524509 


35243859 


17 


121759470 


106779600 


114089288 


105117832 


18 


331072990 


265198150 


304434061 


279216083 


19 


937563530 


840669610 


894584372 


825140032 


20 


2642381430 


2092703550 


2399386239 


2199738033 


21 


7227151280 


6573888100 


6988332717 


6459329037 


22 


20931973090 


16491425740 


18844561759 


17267339059 


23 


55793302330 


51185968460 


54473434666 


50419312152 


24 


164764171030 


129673789110 


147471723662 


135162732506 



Table 2: The number of fixed origin n-step SAWs for various starting points x_l in the 
rhombic Penrose tiling, with starting point coordinates {x,y,z) given in the internal space. 
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3.2 Mean number of walks 

Given that a pair of vertices from A{W) are adjacent if and only if they are adjacent in L, 
we see that the neighbours of a vertex with image x_l can be found by sequentially adding 
the E± image of all possible edges in L to and testing if the new points lie in W. If it 
does the adjacent vertex exists in A{W). By recursively checking all possible neighbours 
of a vertex, all possible walks on the lattices will be found. 

Given an origin vertex in A{W), we know its Ej_ image x^j_ must lie somewhere in 
W, i.e. x*} G W'^ = W {W" is the region x^j_ can lie in given our knowledge of the n steps 
in the walk). If we take a step s (with projection s± onto E±) to a possible adjacent vertex 
x^, then we know x]_ = x']_ + s±. Furthermore if x^ G A(l^) is true, x]_ G W. Hence 
W'^ = (W n + s±)) — s±. The probability that the step s is possible from a random 
x^ is given by the ratio between the areas of and W. 

Extending this to a walk of length n with steps s\i = 1 . . .n, = {W fl {W''~^ + 
Si=i ^±))~ Yli=i the probability of the walk existing is the ratio of the areas of W"- 
and W. For example, consider the shaded in Figure El for the particular walk in the 
Ammann-Beenker tiling which steps west, south, south-west, north-west, west, north then 
north. For the Ammann-Beenker tiling, these probabilities are of the form a + bX, where 
A = 1 + ^2 and a, 6 G Q. 

Adding the self-avoiding constraint and summing the probabilities results in the ex- 
pected number of SAWs beginning at a random origin. Again note that applying this 
method to a regular lattice would result in the usual SAW series. Due to the extra com- 
plexity of the rhombic Penrose tiling acceptance window, the area calculations are more 
involved, see also [2H]- They lead to mean numbers of the form a + br, where r = (l + -\/5)/2 
is the golden number, and a,b & "Z. The rhombic Penrose tiling also allows steps in ten 
directions, more than the Ammann-Beenker tiling's eight. These facts combine to allow 
greater length series to be computed on the Ammann-Beenker tiling. 
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Figure 5: Examples of . . . for a particular walk on the Ammann-Beenker tiling. 



n 


Ammann-Beenker 


rhombic Penrose 





1 


1 


1 


4 


4 


2 


52-16A 


62 -30r 


3 


80-16A 


-4 + 28r 


4 


444-134A 


914 -488t 


5 


1280-380A 


-820 + 732r 


6 


4492- 1430 A 


13842 -7894r 


7 


10848-3248A 


-17732 + 12860r 


8 


60988-21700A 


173876 -101988r 


9 


89800-27036A 


-255784 + 173720r 


10 


643248-237732A 


1923078 -1143988r 


11 


979776-324200A 


-3149856 + 2073192r 


12 


5486960-2043420A 


19566548 -11734706r 


13 


10785736-3819788A 


-34951044 + 22612992r 


14 


45253532-16927618A 


192557132 -116151274r 


15 


110294592-40576780A 


-366912524 + 234803904r 


16 


375796808-141368464A 




17 


1058437232-398339560A 




18 


3259350860-1238175678A 




19 


9526156024-3632872284A 




20 


29127575440-11192322668A 




21 


81536068712-31337365980A 




22 


259724099656-100797073134A 





Table 3: The mean number of n-step SAWs for the Ammann-Beenker tiling and the rhombic 
Penrose tiling, where A = 1 -|- \/2 and r = (1 + x/5)/2. 
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3.3 Total number of walks 



Investigating all possible walks as in the mean number of walks method, we count instead 
the number of non-zero contributions to the mean value. This counts the number of trans- 
lationally inequivalent walks with having positive area or, equivalently, the number of 
translationally inequivalent walks which may occur anywhere in the tiling. Again we note 
that applying this method to a regular lattice gives the usual SAW series. 



n 


Ammann- Beenker 


Penrose 





1 


1 


1 


8 


10 


2 


56 


90 


3 


288 


560 


4 


1280 


2800 


5 


5344 


12060 


6 


20288 


48520 


7 


74192 


182000 


8 


260336 


658300 


9 


892800 


2282400 


10 


2976512 


7749440 


11 


9828256 


25634920 


12 


31758112 


83615140 


13 


101847216 


268113660 


14 


322240144 


850895040 


15 


1012048208 


2668534600 


16 


3147031584 




17 


9732815728 




18 


29852932384 




19 


91182029360 




20 


276695822928 




21 


836719766336 




22 


2516664888416 





Table 4: The total number of n-step SAWs for the Ammann-Beenker tiling and the rhombic 
Penrose tiling. 



3.4 Self-avoiding polygons 

A self-avoiding polygon is equivalent to a self-avoiding walk in which the initial and final 
vertices are adjacent. In the enumeration of self-avoiding polygons, we do not distinguish 
between polygons having, up to a translation, the same boundary but different fillings of 
the interior. 
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SAPs may be enumerated in the same manner as we enumerate SAWs. The additional 
property of end point adjacency allows the backtracking algorithm to be pruned earlier. 
When enumerating SAPs up to size N, a walk that visits a vertex after D steps that is 
further than N — D steps from the origin may be pruned from the search tree. Such a 
walk can never form part of a SAP of < steps. This was used to extend the length of 
the rhombic Penrose SAP series. The Ammann-Beenker SAP series were calculated at the 
same time as the SAW series and hence are of the same length. The extensive run time 
requirements precluded further series extension. 

For the computation of occurrence probabilities of self-avoiding polygons, the loop 
vertices are taken into account, as described in Section 3.2 for the mean number of walks. If 
the self-avoiding polygon has n loop vertices G A{W), where i = 1, . . . ,n, the acceptance 
domain is W"' = C\ii^ ~^±)y occurrence probability is given by the ratio between 

the areas of W"" and W, see also |2H]- 



n 


mean number 


total number 


mean number 


total number 


2 


4 


8 


4 


10 


4 


8 


48 


8 


80 


6 


12A 


384 


108-48r 


840 


8 


800-272A 


2960 


240-64r 


6480 


10 


2840-880A 


21600 


6192-3364r 


49760 


12 


28152-9984A 


170256 


25584-13248r 


394080 


14 


47712-9884A 


1322048 


179200-95340r 


3087140 


16 


869600-299392A 


10194720 


162976-5440r 


24020160 


18 


215712+294408A 


79960896 


2704140-1067580r 


183529440 


20 


14980920-3730840A 


618248240 






22 


152588920-47048100A 


4726263168 







Table 5: The mean number of ra-step SAPs and the total number of SAPs for the Ammann- 
Beenker tiling (first two columns) and for the rhombic Penrose tiling (last two columns), 
where A = 1 + ^2 and r = (1 + V5)/2. 



4 Analysis of series 

The various sequences were analysed using standard methods of asymptotic analysis of 
power series expansions as described in ^j. For self-avoiding walks and polygons, it is 
easy to prove that the limit lim„^oo (cn)^''" exists by use of concatenation arguments PE] . 
We assume the usual asymptotic growth of the sequence coefficients c^, viz: 

Cn = Ax-'^n^-^ [I + 0{n~')\ (n ^ oo, < e < 1). (4) 

On the square lattice, there is overwhelming evidence jH] of the above asymptotic behaviour 
with 7 = 43/32. There is however no proof of this assumption. For interesting new 
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Figure 6: Ratio plot of Cn/cn-i against 1/n for fixed origin (full circles) and mean number 
(large empty circles) Ammann-Beenker SAWs. The square lattice SAW data is plotted in 
small circles for comparison. 



developments see j2S] • The above assumption results in an asymptotic growth of the ratios 



Cn 1 
•-n - - — 



1 + ^-^ + 0{n~'~') 
n 



(n ^ 00,0 < e < 1), (5) 

Cm— 1 -^c 

which may be used to extrapolate numerical estimates of Xc and 7. Whereas it has been 
proved for the square lattice that the limit lim„^oo Cn/c„_2 exists and coincides with x~'^ 
, a similar statement for the ratios r„ is not known. For some lattices, counterexamples 
are known [13j. 

Fig. ini shows a plot of the ratios r„ against 1/n for a typical fixed origin Ammann- 
Beenker walk (full circles) and for the ratios of the mean numbers of Ammann-Beenker 
walks (large empty circles). We notice that the fixed origin data suffers from dramatic fluc- 
tuations, which are smoothed out by averaging, but are still larger than the corresponding 
square lattice data |H|, which is shown in small circles. The oscillating behaviour of the 
mean number of walks data is due to an additional singularity of the sequence generating 
function at a; = —Xc, which, for the case of the square lattice, is well understood due to 
anti-ferromagnetic ordering [H]. To obtain estimates of Xc and 7, we used the standard 
method described in p2] and first mapped away the singularity on the negative real axis 
by an Euler transform and then used Neville-Aitken series extrapolation. 

We also used the method of differential approximants (DAs) ^2]- The underlying 
idea is to fit a linear differential equation with polynomial coefficients to the generating 
function of the sequence, truncated at some order hq. The critical points and critical 
exponents of the differential equation are expected to approximate the critical behaviour 
of the underlying sequence. Application of first order and second order DAs has proved 
useful for the analysis of square lattice SAWs [T^. A first order DA involves fitting the 
coefficients to the differential equation Pi{x)xf'{x) + Po{x) f {x) = R{x) where Pi{x), Po{x) 
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and R{x) are polynomials of degree N, M, and L respectively. We refer to this as a 
[L/M;N] DA. We analysed the approximants [L/N -l;N], [L/N;N], [L/N + I; N] for 
1 < L < 8. We computed estimates for Xc and 7 by averaging of the different DA results, 
given a fixed number of series coefficients hq. Since this yields more accurate estimates 
of Xc and 7 (typically of one digit better) than the ratio method and series extrapolation, 
we list in Table IHl only the results for the DA analysis. Note that the errors are no strict 
error bounds, but arise from averaging over approximants [L/Nq]Ni] for different values 
of A'^o, A^i and L, as explained in [12j. 





mean no. 


(0,0) 


/I V2 V2\ 


(1,0) 


(1- v^,0) 




0.36414(18) 


0.3659(14) 


0.3644(13) 


0.3644(21) 


0.3657(16) 


7 


1.325(19) 


1.45(16) 


1.30(14) 


1.32(17) 


1.46(17) 



Table 6: Estimates of Xc and 7 for Ammann-Beenker SAWs. Numbers in brackets denote 
the uncertainty in the last two digits. 

As suggested by the ratio plot, the estimate using the data for the mean number of 
walks yields the most precise estimates, which are, however, one order of magnitude in 
error worse than the corresponding estimates for square lattice SAWs. 

An analysis of the total number of SAWs on the Ammann-Beenker tiling using first 
order DAs yields Xc = 0.3647(33) and 7 = 3.14(37). Whereas the critical point estimate 
is consistent with the previous analysis, the exponent estimate deviates from the value of 
43/32 = 1.34375 for fixed origin SAWs or the mean number of SAWs. This phenomenon 
reflects the fact that the number of Ammann-Beenker patches of radius r grows asymp- 
totically as r^. (More generally, for aperiodic Delone sets in M.'^ described by a primitive 
substitution matrix, the number N{r) of patches of radius r grows like N{r) ~ r'^ fl^\'21\.) 
Since the SAW has fractal dimension 4/3, we expect an asymptotic increase of the number 
of SAWs by n'^", where u = 3/4. Thus 7 = 43/32 + 2z/ = 2.84375. Data extrapolation is 
consistent with this value. 

The analysis of SAP data follows the same lines. However, the estimates suffer from 
large finite size errors due to the low number (11) of available coefficients. First order 
differential approximants for the mean number of SAPs yield Xc = 0.3688(41). We assume 
Xc{SAP) = Xc{SAW), which has been proven for the square lattice case |14j. 

For the critical exponent a = 2 + 7, we expect for the mean number of SAPs by 
universality that a = 1/2, being the believed exact value for the square lattice (and 
numerically confirmed to very high precision ^IBI). Due to lack of data it is not possible 
to give estimates of critical exponents. An analysis of the total number of SAPs on the 
Ammann-Beenker tiling using first order DAs yields Xc = 0.3587(15). 

The above analysis has also been applied to the rhombic Penrose tiling data. We observe 
qualitatively the same finite size behaviour as for the Ammann-Beenker tiling data, though 
the fluctuations are a bit less pronounced. In Table[7|we list estimates for Xc and 7 obtained 
by analysing first order differential approximants. 

An analysis of the total number of SAWs on the rhombic Penrose tiling using first 
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mean no. 


(0,0,0) 


(0,0,1) 


(0.25,0.5,0) 


(0.5,0.5,1) 




0.36322(29) 


0.3621(12) 


0.3613(16) 


0.36248(83) 


0.36347(51) 


7 


1.333(26) 


1.28(14) 


1.19(22) 


1.303(83) 


1.387(62) 



Table 7: Estimates of Xc and 7 for Penrose SAWs. Numbers in brackets denote the 
uncertainty in the last two digits. 

order DAs yields Xc = 0.3638(31) and 7 = 2.77(21). The estimate of the critical point is 
consistent with the estimates from the other methods of counting. For the critical exponent, 
we again expect a value of 7 = 43/32 + 2z/ = 2.84375, which agrees with the extrapolation 
within numerical accuracy. 

For the analysis of SAP data, only 9 series coefficients are available. First order dif- 
ferential approximants for the mean number of SAPs yield Xc = 0.372(11). An analysis 
of the total number of SAPs on the rhombic Penrose tiling using first order DAs yields 
Xc = 0.3590(22). Again, due to lack of data, it is not possible to extrapolate reasonable 
estimates for the critical exponent a. 

Conclusion 

We extended previous enumerations for self-avoiding walks and polygons on the Ammann- 
Beenker tiling and on the rhombic Penrose tiling and extracted estimates for the critical 
point and critical exponent, using different counting schemes. It turned out that averaging 
with respect to the occurrence probability in the whole tiling leads to the best estimates, 
whereas data produced by fixing an origin leads to strong finite size oscillations. The results 
support the universality hypothesis that the critical exponents appear to be the same as 
for the square lattice, within confidence limits. For the total number of walks (polygons) 
we obtain a new exponent refiecting the polynomial complexity of the number of patches 
of the underlying tiling. 

Since the results were obtained using the enumeration method of backtracking, one 
might ask if more efficient enumeration methods can be applied in order to substantially 
increase the length of the series, and hence the accuracy of the estimates. Unfortunately, 
the successful finite-lattice method cannot be applied in this case, without substantial 
development. 

On a mathematically rigorous level, it may be possible to show the equality of the 
critical points for SAPs and SAWs on quasiperiodic tilings by appropriately modifying the 
existing proofs for the hypercubic lattice ^3]. Furthermore, it would be interesting to 
carry out an analysis to determine if random walk behaviour can be proved for dimensions 
greater than four |26j . 

Self-avoiding polygons may also be counted by area. Since this leads to a three-variable 
generating function due to the different areas of the prototiles, it is tempting to ask whether 
the scaling behaviour of these objects is different from that recently found [^01 for self- 
avoiding polygons on two-dimensional lattices. 
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